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Abstract 

Reporting about the formalism with the Dirac equation we describe the dynamics of chiral 
oscillations for a fermionic particle non-minimally coupling with an external magnetic field. For 
massive particles, the chirality and helicity quantum numbers represent different physical quantities 
of representative importance in the study of chiral interactions, in particular, in the context of 
neutrino physics. After solving the interacting Hamiltonian (Dirac) equation for the corresponding 
fermionic Dirac- type particle (neutrino) and quantifying chiral oscillations in the Dirac wave packet 
framework, we avail the possibility of determining realistic neutrino chirality conversion rates by 
means of (helicity) polarization measurements. We notice that it can become feasible for some 
particular magnetic field configurations with large values of B orthogonal to the direction of the 
propagating particle. 
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In spite of their electric charge neutrality, neutrinos can interact with photons through 
loop (radiative) diagrams which represent the interaction between a fermionic field ip(x) and 
an electromagnetic field A^ix). The correspondent Lagrangian term is written in terms of 
the field-strength tensor F^{x) = d»A v {x) - d u A»(x) as 

£ = I ^ F ^(x) - df u (x)\ ^(x) + h.c. (1) 

where x = (t,x), a^ u = f [7^,71/], the dual field-strength tensor T^ v (x) is given by T >xv {x) = 
|e^ A<5 F x5 (x)[3A] and the coefficients /j, and d respectively represent the magnetic and the 
electric dipole moment which establish the neutrino electromagnetic coupling. Circumvented 
by this extensively explored scenario [1-5], our aim is to partially accommodate the Dirac 
wave packet formalism [6, 7] with this further class of static characteristics of neutrinos, 
namely, the (electro)magnetic moment which appears in a Lagrangian with non-minimal 
coupling. By following this line of reasoning, we evaluate the corresponding Hamiltonian 
equation, from which we can describe the dynamics of the chiral projection operator (75) 
and, consequently, determine the complete expression for the neutrino chirality conversion 
rate. Thus we analyze the quantum oscillation phenomena for some particular magnetic field 
configurations with large values of B parallel to the direction of the propagating particle. 
We firstly set the non-minimally interacting (Dirac) Hamiltonian 

H = a-p + 0[m- (i(S- B(x) + dY.- E(x)], (2) 

in terms of the Dirac matrices, a = Yll=i a ^ — Yll=i lolkk, (3 = 7 , where B(x) and E{x) 
are respectively the magnetic and electric fields. In fact, the Eq. (2) could be extended 
to an equivalent matrix representation with flavor and mass mixing elements {m hj ) where 
the diagonal (off-diagonal) elements described by fi^ and d itj , where i, j are mass indices, 
would be called diagonal (transition) moments. In addition, for both Dirac and Majorana 
neutrinos, we could have transition amplitudes with non-vanishing magnetic and electric 
dipole moments [8-10]. Otherwise, the CP invariance holds the diagonal electric dipole 
moments null [11]. Particularly for Majorana neutrinos, it can be demonstrated that the 
diagonal magnetic and electric dipole moments vanish if CPT invariance is assumed [8]. 
Turning back to the simplifying example of diagonal moments and assuming CP and CPT 
invariance, we can restrict our analysis to the coupling with only an external magnetic field 
B(x) by setting d — 0. From this point, the expression for n can be found from Feynman 
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diagrams for magnetic moment corrections [11] and turns out to be proportional to the 
neutrino mass (matrix), /i = g ^^ 2 m = 2.7 x 1CT 10 fi B ^ where G is the Fermi constant and 
ttijsj is the nucleon mass [35]. 

Assuming the simplifying hypothesis of a uniform magnetic field B, the physical impli- 
cations of the non-minimal coupling with an external magnetic field can then be studied by 
means of the eigenvalue problem expressed by the Hamiltonian equation 

H(p) <p n = E n (p) <p n = {a ■ p + [m - n S • B]} ip n . (3) 

The most general eigenvalue (E n (p)) solution of the above problem is given by 



E n (p) = ±E s (p) = ±yj 'm 2 + p 2 + a 2 + (-iy2^m 2 a 2 + (p x^, s = 1, 2, (4) 

where we have denoted £ , n=1)2)3) 4 = ±^ s=lj2 , and we have set a = fiB. The complete 
set of orthonormal eigenstates tp n thus can be written in terms of the eigenfunctions U(p s ) 
with positive energy eigenvalues (+E 3 (p)) and the eigenfunctions V(p 3 ) with negative energy 
eigenvalues (—E s (p)), 

\ t 



A- a+ A-oti 





t 



^)--^){V^-#'-v^'" 1 )' (5) 

where p s is the relativistic quadrimomentum, p 3 = (E s (p),p), N(p s ) is the normalization 
constant and Af = A 2 (p) ± 2m\a\ - a 2 , af = 2E s {p)\a\ ± (A 2 (p) + a 2 ) with A 2 (p) = 
E 2 s {p) — (m 2 + p 2 ) + a 2 . We can observe that the above spinorial solutions are free of any 
additional constraint, namely, at a given time t, they are independent functions of p and they 
do not represent chirality/helicity eigenstates [15-17]. The time-evolution of a plane-wave 
packet ip(t,x) can be written as 

i;(t,x) = [J£-J2{b(jP.)U(jp t ) exp[-iE s (p)t]+d*(p s )V(p s ) exp [+iE s (p) t] } exp [i p ■ x] , (6) 



with p s = (E s ,—p). Meanwhile, the Eq. (6) requires some extensive mathematical manip- 
ulations for explicitly constructing the dynamics of a generic operator 0{t). If, however, 
the quoted observables like the chirality j 5 , the helicity h or even the spin projection onto 
B commuted with the Hamiltonian H, we could reconfigure the above solutions to simpler 
ones. Let us then assume that the magnetic field B is either orthogonal or parallel to the 



momentum p. For both of these cases the spinor eigenstates can then be decomposed into 
orthonormal bi-spinors as y?* 2 and xt, 2 which are eigenstates of the spin projection operator 
cr ■ B, i. e. beside of being energy eigenstates, the general solutions U(p 3 ) and V(p 3 ) would 
become eigenstates of the operator S ■ B and, equivalently, of S • a. Now the Eq.(3) can be 
decomposed into a pair of coupled equations like 



(±E a -m + a - a) <pf 
(±E, +m-a-a) X ^ 



±<x -Pxt, 
±<r-P<pf, 



(7) 



where we have suppressed the p 3 dependence. 

By introducing the anti-commuting relation {er ■ p, cr ■ B} = derived when p ■ B = 0, 
the eigenspinor representation can be reduced to 

" cr p 



W(p.) 



e + m 
2e 



cr p + 

eo+m s 



and V(p s ) = 



e + m 
2e 



eo+m 



x; 
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with s = ^/p 2 + m 2 and the energy eigenvalues ±E S = ± [e + (-l) s |a|]. 

Analogously, by introducing the commuting relation [cr • p, cr • B] = when p x B = 0, 
the eigenspinor representation can be reduced to 



2E, 



" p 

E s +m s rs 



and V(p s ) 



E. 



2E S 



a p 

Es+m s 



(9) 



with m 3 = m — (-l) s |a| and the energy eigenvalues ±E 3 = i^/p 2 + m 2 s . Since we can set 
ip+ 2 = Xi 2 for both of the above cases as the components of an orthonormal basis, the 
orthogonality and the closure relations are immediate. 

Finally, the calculation of the expectation value of 7 5 (i) is substantially simplified when 
we substitute the above closure relations into the wave-packet expression of Eq. (6). To 
clarify this point, we suppose the initial condition over ip(t,x) can be set in terms of the 
Fourier transform of the weight function (p(p - p ) w = XL =i 2 {b(Ps)U(p s ) + d*(ps)V(p s )}, 

dp 



(27T) 



r^O^Po) ex P [ip-x]w 
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where w is some fixed normalized spinor, and we have found [6] b(p 3 ) = tp(p-p a )W{p s )w 
and d*(p s ) = ip(p~ p )V j; (p a )w. For any initial state ip(Q,x) given by Eq. (10), the negative 
frequency solution coefficient d*(p s ) necessarily provides a non-null contribution to the time- 
evolving wave-packet. This obliges us to take the complete set of Dirac equation solutions to 



construct a complete and correct wave-packet solution. Only if we consider the initial spinor 
w being a positive energy (E s (p)) and momentum p eigenstate, the contribution due to the 
negative frequency solutions d*(p s ) will become null and we will have a simple expression 
for the time-evolution of any physical observable. By substituting the explicit form of the 
spinors of Eqs. (8) and (9) (closure relations) into the time-evolution expression for the 
above wave-packet, the Eq.(6) can thus be rewritten as 

W,x)= [J^-pip-pJexpiip-x] {[cos[£?.t]-if sin [£,*]] ( M-D^oS-A ) j w (n) 

^ ' 3 = 1,2 

for the first case where E s = [s + (-l) s |a|] and H = ct ■ p + 7 m, or as 

W,x) = J-^<p(p-p )exp[ip-x] {[cos[E s t]-i^sm[E s t]\ ( lzL=2^A "j j w (12 ) 

for the second case where i? s = y'p 2 + m 2 and H s = a. ■ p + 7 7?v 

It was already demonstrated that, in vacuum, chiral oscillations can introduce very small 
modifications to the neutrino conversion formula [21, 22]. The probability of a neutrino 
produced as a negative chiral eigenstate be detected after a time t can be summarized by 
the expression 

P{y a ,L ^v a ,L\t) = j ' d>xi;Ht,x)^^i>(t,x) = ^(l-( l5 )(t)). (13) 

From this integral, it is readily seen that an initial -1 chiral mass-eigenstate will evolve with 
time changing its chirality. By assuming the fermionic particle is created at time t = as a 
-1 chiral eigenstate (^ 5 w = -w), in case of {cr -p, a ■ B} — (B orthogonal to p), we could 
write 

<75>(*) = (-1) /(|^ 2 (P-Po){f cos[2|a|i] + ^cos[2e t]}, (14) 
where we have used the wave packet expression of Eq. (11) and, in addition to w^ 5 w — -1, 
we have also observed that {H , 7 5 } = 27 5 S ■ p and, subsequently, w^H • P7 ( ,S ■ aw = 0. 

2 

For the above result, in addition to the non-interacting oscillating term cos [2 e £], which 
comes from the interference between positive and negative frequency solutions of the Dirac 
equation, we have an extra term which comes from the interference between equal sign 
frequencies and, for very large time scales, can substantially change the oscillating results. 
In this case, it is interesting to observe that the ultra-relativistic limit of Eq.(14) leads to 
the following expressions for the chiral conversion formulas, 

P{v a ,L ~ \ (! + (-) cos [2 \a\t\) (15) 



which, differently from chiral oscillations in vacuum, can be phenomeno logically relevant. 
Obviously, we are reproducing the consolidated results already attributed to neutrino spin- 
flipping [11] where, by taking the ultra- relativistic limit, the chirality quantum number can 
be approximated by the helicity quantum number, but now it was accurately derived from 
the complete formalism with Dirac spinors. We still remark that, in the standard treatment 
of vacuum neutrino oscillations, the use of scalar mass-eigenstate wave packets made up 
exclusively of positive frequency plane-wave solutions is implicitly assumed. 

Turning back to the case where [a ■ p, <r ■ B\ = (B parallel to p), we could have 
a phenomeno logically more interesting result. By following a similar procedure with the 
mathematical manipulations, we could write 

= /70i^(p-^)E{[S + i cos [ 2 ^^ ^( M -T s ' d H, (i6) 

^ ' 3 = 1,2 

where we have used the wave packet expression of Eq. (12) and, in the second passage, 
we have observed that w^ 5 w = -1, w^[H s , ^ 5 ]w = and H s ^ r> H a = p 2 — m 2 . The above 
expression can be reduced to a simpler one in the non-interacting case where just the mass 
term appears. Due to a residual interaction with the external magnetic field B, we could 
also observe chiral oscillations in the ultra-relativistic limit. As a toy model illustration, by 
assuming a highly peaked momentum distribution centered around a non-relativistic mo- 
mentum Pi <C m 3 , where the wave packet effects are practically ignored, the chiral conversion 
formula can be written as 

P{v a ,L —> ^a,R',t) ~ i ^1 — cos [2 mi] cos [2 \a\ t] — sin [2 mi] sin [2 \a\ t]w^H ■ av?j (17) 

where all the oscillating terms come from the interference between positive and negative 
frequency solutions which compose the wave packets. 

By assuming that the neutrinos before interacting with the magnetic fields are produced 
with completely random spin orientation, which can be viewed as a collection of neutrinos 
where 50% are characterized by spin-up states and the remaining 50% by spin-down states, 
as we have noted before, we can understand the result provided by the Eq. (16) as a phe- 
nomenological instrument much more interesting than the result expressed in the Eq. (14). 
By averaging on time both the above results obtained for (75) (t), we can easily observe that 
the corresponding average value for the case where p is perpendicular to B is zero, and the 
same average value for the case where p is parallel to B depends on the polarization of the 
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initial state as we can observe in the following expression, 

«7 5 >(*)We = (-i) [^ 2 (p-po) £ [% w] ( M -T S '* H (18) 

^ ' 3 = 1,2 

By supposing that the construction of a detection apparatus which identify the neutrino 
state of polarization (helicity!) is feasible, the last result becomes phenomenologically rel- 
evant since it distinguish the chirality conversion rates for spin-up (s = 1) and spin-down 
(s = 2) states. Once we have assumed the neutrino electroweak interactions at the source 
and detector are {left or negative) chiral (^^(l — ^ipW^) only the component with neg- 
ative chirality contributes to the effective detection result. The residual contribution due 
to the chiral oscillation mechanism appears in the final time averaged result. It leads to 
different expected probabilities of detecting a negative chiral eigenstate for the opposite po- 
larization states. By comparing the initial polarization states where, before interacting with 
the magnetic fields, for any detection process, the spin-up and spin-down states would be 
equally probable, with the expected results for the polarization states after traversing the 
magnetized region with p parallel to B, where the chiral oscillation mechanism effectively 
takes place, we observe that the chiral oscillation mechanism can lead to a modification 
on the final expectation values. It reflects in a ratio between the spin-up and spin-down 
negative chiral eigenstate probabilities different from the unity as we have illustrated in the 
Figs. 1 and 2. 

2 

As one can notice in the plots, since the coefficient of the oscillating term goes with , 
the chiral oscillation effects are not relevant for ultra-relativistic neutrinos [22, 23]. At the 
same time, there is a natural scale for the magnetic field strength that is required to have 
a significant impact on the quantum process, i. e. a critical value of 4.41 x 10 13 G. There 
are reasons to expect that magnetic fields of such or even larger magnitudes can arise in 
cataclysmic astrophysical events of such as supernova explosions or coalescing neutron stars, 
situations where an enormous neutrino outflow should be expected. Two classes of stars, 
the soft gamma-ray repeaters (SGR) [27] and the anomalous x-ray pulsars (AXP) [28, 29] 
are supposed to be remnants of such events which form magnetars [26], neutron stars with 
magnetic fields of 10 14 — 10 15 G. From the theoretical point of view, the possibility of larger 
magnetic fields 10 16 — 10 17 G have not been discarded yet. The early universe between the 
QCD phase transition (~ 1CT 5 s) and the nucleosynthesis epoch (10~ 2 — 10 2 s) is believed to 
be yet another natural environment where strong magnetic fields and large neutrino densities 
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FIG. 1: We plot the ratio between the detection probabilities of spin-up and spin-down nega- 
tive chiral eigenstates as a function of the magnetic interaction energy /j,B/m for different values 



with a completely random spin orientation, which can be viewed as a collection of neutrinos where 
one-half are characterized by spin-up states and the remaining one-half by spin-down states. 

could coexist. 

In which concern with the polarization measurements, the observation of the dependence 
of the neutrino chirality conversion rate on the neutrino polarization state (measurements) 
can be converted into a clear signal of the presence of right-handed (positive chiral) neutrinos 
in the neutrino-electron scattering. In fact, in a more extended scenario, the scattering of 
neutrinos on a polarized electron target [24] was proposed as a test for new physics beyond 
the Standard Model (SM). To search for exotic right-handed weak interactions [24, 25] the 
strong polarized neutrino beam and the polarized neutrino target is required [24, 25]. It 
has been shown how the presence of the right-handed neutrinos changes the spectrum of 
recoil electrons in relation to the expected standard model prediction, using the current 
limits on the non-standard couplings. In this framework, the interference terms between the 
standard and exotic couplings in the differential cross section depend on the angle between 
the transverse incoming neutrino polarization and the transverse electron polarization of the 
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FIG. 2: We plot the ratio between the detection probabilities of spin-up and spin-down negative 
chiral eigenstates as a function of p /m representing the propagation regime for different values 

of the magnetic interaction energy fiB/m. We have also established that ^ - — ^ — ^ u; = 

+ /i-(-i) 2 £-d\ i 
™ T J w = 2- 

target, and do not vanish in the limit of massless neutrino. 

Finally, we would have been dishonest if we had ignored the complete analysis of the 
general case comprised by Eqs. (3-5) where we had not yet assumed an arbitrary (simplified) 
spatial configuration for the magnetic field. We curiously notice the fact that those complete 
(general) expressions for propagating wave packets do not satisfy the standard dispersion 
relations like E 2 = m 2 + p 2 excepting by the two particular cases where E s (p) 2 = m 2 s +p 2 for 
p x B = 01 e' 2 = m 2 + p 2 for p ■ B = 0. In particular, the process of neutrino propagation 
through an active medium consisting of magnetic field and plasma and the consequent 
modifications to the neutrino dispersion relations have been studied in the literature [30]. 
In addition, such a general case leads to the formal connection between quantum oscillation 
phenomena and a very different field. By principle, it could represent an inconvenient 
obstacle forbidding the extension of these restrictive cases to a general one. However, we 
believe that it can also represent a starting point in discussing dispersion relations which 
may be incorporated into frameworks encoding the breakdown (or the validity) of Lorentz 
invariance. 
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Turning back to the foundations for applying the intermediate wave packet formalism in 
the neutrino oscillation problem, we know the necessity of a more sophisticated approach is 
required. It can involve a field-theoretical treatment. Derivations of the oscillation formula 
resorting to field-theoretical methods are not very popular. They are thought to be very 
complicated and the existing quantum field computations of the oscillation formula do not 
agree in all respects [31]. The Blasone and Vitiello (BV) model [32, 33] to neutrino/particle 
mixing and oscillations seems to be the most distinguished trying to this aim. They have 
attempt to define a Fock space of weak eigenstates to derive a nonperturbative oscillation 
formula. Also with Dirac wave-packets, the flavor conversion formula can be reproduced [7] 
with the same mathematical structure as those obtained in the BV model [32, 33]. 

To summarize, in the context where we have intended to explore the Dirac formalism, 
we have assumed that a satisfactory description for understanding chiral oscillations of 
fermionic (spin one-half) particles like neutrinos requires the use of the Dirac equation as 
evolution equation for the mass-eigenstates. With such a framework stated, we have added 
the non-minimal coupling of the neutrino magnetic moment with an electromagnetic field in 
an external interacting process [1]. Although clear experimental evidences are still missing, 
we have reinforced the idea that the spinorial form and the interference between positive 
and negative frequency components of the mass-eigenstate wave packets can introduce small 
modifications to the standard conversion formulas where chiral oscillations are taken into ac- 
count, in particular, under circumstances where we treat the phenomeno logical possibility 
of neutrino interactions with very large magnetic fields. For the future outline, we intend 
to inspect the correlated constructions with more general magnetic field spatial configura- 
tions as well as we plan to aggregate the corrections discussed here to a more generalized 
study of chiral conversion rates, polarization effects and neutrino propagation in magnetized 
media which have been already under consideration in the literature, in particular, those 
which could be interesting from the point of view of observable effects or phenomenological 
implications. Also in the theoretical framework, we can notice a more extended discussion 
involving the evolution equation of the covariant neutrino spin operator in the Heisenberg 
representation in presence of general external fields [5]. Such an analysis involves both elec- 
tromagnetic and weak interactions with background matter [4, 5]. In a similar context, some 
developments considering the minimal coupling of an electrically charged particle with an 
external magnetic field were also performed for an electron described by the Dirac equation 
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[3] and the theory of spin light of neutrino in matter and electromagnetic fields has been 
extensively studied [5]. 
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